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Abstract 


A new approach is developed to evaluate contribution of the electron-electron correlations into 
bremsstrahlung from few-electron ions and atoms. Our approach is based on the explicit formula 
for the electron density distribution in such systems. We derive the closed analytical formula for 
the matrix elements which are needed for highly accurate computations of atomic form-factors of 
two-electron atoms and ions. We also discuss the energy loss due to bremsstrahlung in a plasma 
which contains multi-charged ions and free electrons. Bremsstrahlung from a high-temperature 
plasma is considered as well as its role in the high-temperature burn-up of deuterium plasma. 
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I. INTRODUCTION 


As is well known electron-electron correlations in atoms and ions contribute to many 
bound state properties of these systems. Recently, it appears that such correlations also 
affect a number of properties of atoms and ions, including the influence of contributions 
from unbound atomic states. In particular, in this study we consider the braking radiation 
(or bremsstrahlung) emitted by a fast electron when it is slowed down, or even stoped, 
by a few electron ion/atom. In general, this process is represented by the equation [l|: 

+ e~ = + e~ + hu, where the notation A^^ stands for a multi-charged (atomic) 

ion which has a positive electric charge q+, while the notations e~ and fioj designate the 
fast electron and emitted photon, respectively. It appears that the bremsstrahlung cross- 
section and corresponding energy loss explicitly and substantially depend upon the electron- 
electron correlations in the incident ions/atom. Our goal in this short study is to investigate 
this effect and its role for accurate evaluations of the cross-section and energy loss due to 
bremsstrahlung in few-electron ions/atoms. Below, we apply the relativistic units in which 
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is the dimensionless fine- 


h = 1, c = 1 and rrie = where a = ^ 
sructure constant. The unit of length in these units coincides with the Compton wavelength 
of the electron Ag, where = aao ^ 3.8615926764-10“^^ cm and uq = 5.2917721067- 
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cm 


is the Bohr radius. 


In this study the braking radiation emitted during a Coulomb collision of a fast 
electron with the positively charged atomic nuclues is designated by the German word 
‘bremsstrahlung’ [l|. Bremsstrahlung has been analyzed with the methods from Classical 
Electrodynamics and later from Quantum Electrodynamics (see, e.g., S. Q and references 
therein). The goal of the QED analysis was to derive the closed analytical formula(s) for 
the bremsstrahlung cross-section as a function of the photon frequency uj and electron’s 
‘quantum numbers’. For a ‘free’ electron a set of ‘good (or ‘convenient’) quantum numbers 
includes the kinetic energy Ek and the vector of momemtum p^, where k = i for the incident 
electron and k = f for the final electron. In actual experiments to measure bremsstrahlung 
intensity it is very hard to measure angular correlations between directions of propagation of 
the incident/final electron and emitted photon. Therefore, for quantum numbers it is better 
to use the kinetic energy of the electron and absolute value of its momemtum Pk =\ Pk \ 
(scalar). Below, we shall assume that we have an electron with the incident quantum num- 
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bers (i?i,pi) which emits one photon hu = u. The hnal quantum numbers of the electron 
are (£' 2 ,P 2 )- In this notation the well known formula for the bremsstrahlung cross-section 
obtained from direct ana^tical QED calculations performed in the Born approximation 
takes the form (see, e.g., [ij], jj]) 
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where m = me is the electron mass, Ae is the electron’s Compton wave length (oq = a“^Ae 
is the Bohr radius) and Q is the electric charge of the atomic nucleus in an ion/atom which 
has Ne bound electrons, where Ne < Q. The factor E[q) used in this formula is the so-called 
form-factor of the ion/atom. This atomic form-factor is dehned traditionally as: 


p-\-oo p p-\-oo 

E{q) = / f ne(r) exp(—q ■ v)r‘^drdVL = dvr / ne{r) sm.{qr)rdr 
Jo J Jo 
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where ne(r) is the electron density in a few-electron atom/ion which is assumed to be 
spherically symmetric. The density ne(r) is normalized and its norm equals the number of 
bound electrons Ne, i.e. ne(r)r‘^dr = Ne- The parameter q in the expression for the 
form-factor E{q) equals to the ratio of the velocity of the fast electron Ve to the atomic 
velocity ac{= 1 in the relativistic units). 

At small non-relativistic energies, when pi <C me = m in relativistic units, the formula, 
Eq.dH), takes the form 
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where Ti = ^ is the non-relativistic kinetic energy of the incident electron. Note that this 
expression does not contain any energy and/or momentum of the hnal electron, i.e any of 
the E 2 and/or p 2 values. In the opposite ultra-relativistic case, when Ei S> m and E 2 3> m, 
the formula for the bremsstrahlung cross-section, Eq. ([!]), is written in the form 
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where fsiq) is the Born anrplitude for elastic scattering of an electron by a Ag-electron 


, e.g., Q and p). In the case of elastic scattering q = 2/i;sin|. This means 


ion/atom (see 

that by determining the scattering amplitude of elastic scattering in the Born approximation 
one also hnds the factor which is needed for calculation of the differential cross-section of 
bremsstrahlung ^ in the ultra-relativistic limit. 

Note that each of the formulas Eq.([T]) and Eqs.([3]) - (jl]) contains the form-factor F{q) 
dehned in Eq.(l2]). For many few-electron ions the overall contribution of the form-factor 
is relatively small and rapidly decreases when the nuclear charge Q grows. Moreover, for 
such atomic systems the electron density ne(r) can be considered (to very high accuracy) 
as a spherically symmetric function. This corresponds to the actual physical picture of 
bremsstrahlung as a braking radiation emitted during Coulomb interaction between a fast 
electron and atomic nucleus. The role of atomic electrons is an electrical screening of the 
central nucleus. In general, atomic form-factors for different atoms/ions can be determined 
to high accuracy by using computational methods of modern atomic physics. Derivation of 
the analytical formulas and numerical calculations of the form-factors for several different 
few-electron ions and atoms is one of the goals of this study. 


II. FORMULAS FOR THE FORM-FACTOR IN FEW-ELECTRON IONS 

As follows from the formulas given in the Introduction to determine the bremsstrahlung 
cross-sections for actual few-electron ions/atoms one needs to evaluate the atomic form- 
factor F{q), Eq.(j2]). In this Section we derive closed analytical expressions for this quantity 
in different few-electron ions/atoms. For one-electron atomic systems we can use hydrogenic 
wave functions. This drastically simplihes all calculations of the atomic form-factor for 
different bound states, i.e. states with the different values of the angular momentum L 
and the principal quantum number n. To simplify the problem let us consider the ground, 
doublet state in hydrogen-like, one-electron ions with nuclear charge Q. In this case 
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the form-factor F{q) is 
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where we used the formula Eq. (3.944.5) from j^. This formula can easily be generalized 
to the case of two-electron ions/atoms, if it is possible to neglect the electron-electron cor¬ 
relations in such systems. In the lowest order approximation such a generalization can be 
achieved with the substitution Q ^ Q — ^ in the formula, Eq.dH]). The atomic form factor 
changes correspondingly. However, the hrst term in the [Q — F{q)]‘^ expression, i.e. the 
electric charge of the nucleus Q, which can be found in each of Eq.([T]) and Eqs.([3]) - (jl]) does 
not change, since Q is the actual nuclear charge in the ions. These arguments lead to the 
following approximate formula for two-electron multi-charged ions 


[Q-F{q)]^ = Q^{l- 
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where A^e = 2 is the total number of bound electrons. These calculations are simple and 
straighforward. However, analogous computations with the use of the truly correlated wave 
functions for few-electron ions/atoms become signihcantly more complicated. 

For relative simplicity, let us consider the truly correlated wave functions of two-electron 
atoms/ions. Such wave functions can be represented in terms of the different variational 
expansions. One of the most effective and accurate expansion is the exponential varia¬ 
tional expansion written in relative/perimetric coordinates. For the ground S'—states this 
expansion takes the form (see, e.g., 8| and S) 

N 

^ = (l + Pu) exp(-Q;ir32 - Arsi - 7ir2i) 

i=l 
N 

= ^Ci exp(-Q(ir32 - Aoi - 7T2i) + exp(-/3ir32 - air3i - 7ir2i) 


( 8 ) 
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The expansion, Eq.([H]), is called the exponential variational expansion in the relative coor¬ 
dinates r 32 ,r 3 i and r 2 i. Each of the three relative coordinates Vij is dehned as the difference 


between the corresponding Cartesian coordinates of the two particles, i.e., Vij =\ fj—r,- |= r 
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It follows from this dehnition that the relative coordinates r 32 , r 3 i and r 2 i are translationally 
and rotationally invariant. Below, the index 3 is used to designate the atomic nucleus, while 
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indexes 1 and 2 stand for atomic electrons. The coefficients Ci are the linear variational 
parameters of the expansion, Eq.dH]), while the parameters ai^jSi and 7 * are the non-linear 
parameters of this expansion. In general, the total energy of the ground S'—state of the 
two-electron ion depends upon the total number of basis functions N used in calculations. 
The operator P 12 in Eq.(|S]) is the permutation operator for two identical particles (elec¬ 
trons). The very high efficiency of the variational expansion, Eq.(| 8 ]), in actual applications 
to the two-electron ions is related to the fact that all non-linear parameters a*, (3i and 7 , are 
carefully varied in such calculations. 

An analytical formula for the form-factor F{q) derived with the use of the exponential 
variational expansion takes the following form 
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First, we can calculate the internal integral in this equation. The result is written in the 
two following forms: ( 1 ) for r 32 > r 3 i 
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I'r32+r3l — r^2 1 \ 

- exp[-(7. + 7,)(r3i 

These two integrals / and J are the functions of the two variables r 32 - 1 - rsi and |r 32 — r 3 i|. 
By performing integration over the r 3 i relative coordinate one obtains the explicit formula 
for the form-factor F{q) as a quantity which depends upon the r = r 32 radial variable (or 
electron-nucleus distance r^N)- In some cases the knowledge of this dependence is crucial to 
understand the nature of the physical problem. 

However, if we need to know only the absolute value of the form-factor F{q), then it 
is possible to use another approach to its calculation. In this case instead of three rela¬ 
tive interparticle coordinates r 32 ,r 3 i and r 2 i we can introduce three perimetric coordinates 
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where Uk = ^{vki + Tkj - r^), Tab = Vba and {i,j,k) = (1, 2, 3) js]. The inverse 
relation take the form r^h = Ua + Ub- The three perimetric coordinates generally are inde¬ 
pendent of each other, non-negative and each of them varies between 0 and -|-cxd. As follows 
from these properties, three perimetric coordinates form a very convenient set of variables in 
order to calculate arbitrary three-particle integrals. Let us apply perimetric coordinates to 
determine the atomic form-factors of the different two-electron ions/atoms. First, consider 
the integral between jo{qr^ 2 ) and i and j exponential basis functions. In relative coordinates 
we can write 
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where the factor jo(?^r 2 ) = is the spherical Bessel function of zero order. In perimetric 

coordinates this integral takes the form 
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sm{qu3) cos(gti2) + cos(gti3) sm{qu2)]{u3U2 + U3U1 + U2U1 + uj) x 


exp[—Z m 3 — Yu 2 — Xui]du 3 du 2 dui 


(13) 


where Z = ai + Uj + (Ji + fJj,Y = ctj -f aj + ^i + and X = A + A + 7* + 7i- factor 2 
in this formula is the Jacobian of the transformation (^ 32 , r 3 i, r 2 i) (ti 3 , M 2 ,'Ui). As follows 
from Eq. flT^ . calculation of the integral I is reduced to the analytical computation of eight 
integrals, where each of the contributing integrals is the product of three one-dimensional 
integrals. For instance, the first contributing integral is 
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The final formula for the (A)-niatrix element of the form-factor F{q) is 
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It is clear that this formula is a regular fuuctiou of q (i.e., it coutaius uo siugularities) aud 
is uumerically stable for arbitrary q. The formula, Eq. fllGp . for matrix elemeuts is similar 
to our formulas derived iu Q for the matrix elemeuts iuvolviug spherical Bessel fuuctious 
which have beeu obtaiued with the use of the same approach. Results of highly accurate 
uumerical calculatious of some atomic form-factors cau be fouud iu Tables I aud II (for 
differeut ious/atoms aud differeut values of q). Table I shows couvergece of the form-factor 
for the grouud RR—state of the °°He-atom determiued with the use of the formula, Eq.(16), 
aud differeut uumber(s) of basis fuuctious N iu Eq.dH]). As follows from Table I couvergeuce 
of the form-factor F{q) upou N iu Eq.(| 8 ]) is very fast for the He atom (grouud S'—state). 
Briefly, we cau say that results from Table I are almost A^—iudepeudeut. The same situatiou 
cau be obesrved with the form-factors of other two-electrou ious showh iu Table II (H“, Ne®"*", 
Ca^®’*', Ni2®+). This iudicates the maiu advautage of our formula, Eq. fllbp . used for uumerical 
calculatious of the form-factors F{q) iu the two-electrou atoms/ious. Formally, it is possible 
to say that our method based ou the formula, Eq. flT 6 p . completely solves the problem of 
highly accurate computatious of the atomic form-factors for two-electrou ious/atoms. 

Table II coutaius form-factors F{q) (iu a.u.) determiued for a uumber of two-electrou 
ious/atoms iu their grouud l^S*—states. It is iuterestiug to observe chauges of the form- 
factors for such atomic systems due to chauges iu the electric charge of the uucleus Q. For 
very compact ious, e.g, for the Ca^®’*', Ni 2 ®+ ious, form-factors chauge very little wheu the 
parameter q varies betweeu 0 aud 11 (see Table III). Ou the other haud, for the weakly-bouud 
H“ iou all substautial chauges of the form-factor are located iu this area of q variatious. The 
same is true for the form-factor of the ueutral He atom. 

The third approach for determiuiug the form-factor F(a) is based ou the followiug ap¬ 
proximate formula kuowu from atomic physics (see, e.g., ll|): 

dTT r+°° siu(q'r) 
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(17) 


where we have assumed that the electrou deusity distributiou is spherically symmetric. At 
large values of q the iutegraud iu the right-haud side of Eq. flTTP is a fast oscillatiug fuuctiou. 
Therefore, the rauge of large values of the qr variable does uot coutribute to the form-factor 
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F{q). For small qr values (gr < gao, where oq is the Bohr radius) we can write 

r+oo , [qrY ^[qrY (qr)^ ^ (gr)® (gr)^° 
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Assuming that the electron density ne(r) is normalized to the number of bound electrons 
Nf. we can write 



where r = r^N is the electron-nucleus distance (scalar coordinate). In the last formula 
we also eliminated the factor dvr which is compensated by the corresponding factors from 
angular parts of the wave functions, or electron density distribution. This means that 
numerical computations of form-factors for few-electrons atoms and ions are reduced to 
accurate calculations of the (r^^) expectation values for k = 1,2,3,.... In reality, for two- 
electron ions/atoms it is possible to evaluate atomic form-factors by using a hnite number 
of terms, e.g., three, or four terms, in Eq. flT^ . since in this case for relatively small g(< 1) 
the series, Eq. dT^ . converges very fast. Some of the (r^^) expectation values determined 
for a few two-electron ions/atoms (in atomic units) can be found in Table III. Accurate 
computations of the atomic form-factors for three- and four-electron atoms/ions can be 
perfomed analogously. Details of such calculations and analogous computations of the form- 
factors for three- and four-electron ions/atoms will be discussed in future studies. 


III. ENERGY LOSS DUE TO BREMSSTRAHLUNG 


The improtance of the derivatives ^ of the bremsstrahlung cross-section dehned by Eq. ([T]) 
and Eqs.dSj) - (jlj) follows from the fact that these values are directly related to the radiative 
energy loss, or, in other words, to the energy loss due to bremsstrahlung. Indeed, let us 
assume that a fast electron with the incident quantum numbers (i7i,pi) moves through 
matter which contains Ni ions/atoms per cm?. Such an electron emits a photon hu = u and 
becomes the hnal electron with quantum numbers (772, P 2 ). The average energy loss due to 
bremsstrahlung per cm of electron path x is given by 


/dEi\ rEi-m (i(j r 1 /■■El-m , da \ , 
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( 20 ) 


where the value ^ is the differential cross-section of bremsstrahlung dehned in the In¬ 
troduction. The radiative cross-section cr^ dehnes the total radiative losses related to 
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bremsstrahlung. It is interesting to note that this formula can be applied to a large number 
of actual systems, e.g., to evaluate radiative loss in hot plasmas and evaluate an additional 
heating in nuclear reactors which contain substantial amounts of fast decaying f3~ isotopes. 
In fact, our original interest in this problem was related to numerical evaluation of the energy 
loss due to bremsstrahlung in a sample containing a large number of /3“ decaying nuclei. 
In an actual fuel rod which has been taken from a working nuclear reactor 1 cm? usually 
contains between 1 • 10® and 5 • 10^^ of decaying nuclei. The exact number depends 
upon the age of the reactor, its type, working hystory and other factors. The total number 
of j3^ decaying nuclei is ~ 15 - 20 smaller than the analogous number of nuclei. 

By using the analytical expression for the ^ derivative, Eq.(IT]), we can determine the 
radiative cross-section cxr in Eg. (1201) . The hnal formula is 
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( 21 ) 


where the function L{x) = —dilog{x+l) = Io~^ = ’^ — diloqjx) = 

Li 2 {—x) is proportional to the well known dilogarithm function dehned by Eq.(27.7.1) in 12 |. 
Note that dehnitions of the dilogarithm function in different books and text books differ from 
each other. In particular, our function L{x) from E a.f|21ll coincides almost exactly with the 
dilogaritm Li 2 {x) dehned in Wolfram Math World 13|, i.e. L{x) = Li 2 {—x). Furthermore, 
for our L{x) = Li 2 {—x) function one Ends: 


TT^ (Inx)^ 




( 22 ) 


At small X we have L{x) = x — ^ ^— ^ + ■ ■ ■■ This power-type expansion and Eq.( 

allows one to determine the asymptotics of the cross-section ar at small and very large 
energies of the incident electron. In particular, at small energies, when Ei m, one 
Ends ar = 1 — . In this case the radiative cross-section does not depend 

explicitly upon the energy of the incident electron Ei. At large energies of the incident 
electron, when Ei ^ m, the analogous formula for multicharged ions is written in the form 
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For neutral atoms one needs to apply a diEerent 
formula, since in this case the bremsstrahlung cross-section converges to the constant limit 
when El ^ 00 and does not include any term which logarithmically diverges at large 
energies. 
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Actual high-temperature plasma consists of multicharged ions and free electrons. There¬ 
fore, a number of other processes also contribute to the energy losses in actual experimental 
systems. The most important of such processes is the non-eleastic electron scattering by 
atoms and positively charged ions. Another process is electron-electron scattering which 
leads to energy loss of the fast electron and accelerations of many ‘secondary’ electrons, or 
electrons. A third process is the Compton scattering of photons by free electrons with 
the emission of different photons. In the first order this later process is strictly prohibited 
in QED, but in the second order it is possible. All these processes which compete with the 
bremsstrahlung have been analyzed in the literature. The energy loss (per unit length of 
the trajectory) of the fast electron during its non-elastic collisions with atoms/ions can be 
evaluated from the following approximate formula [f^ : 

-2.a^alNQ^n(^) (23) 


where 71 = is the electron’s 7 —factor, N is the number of ions/atoms per cm? and I ~ 
13.605 eV is the ‘atomic’ ionization potential. The ratio R of the — and 

derivatives is written in the form: R = where A(^ 1600) is a numerical constant which 
is chosen to produce the best fits for known experimental data. 

Consider now the process of electron-electron scattering which leads to the formation of 
the secondary accelerated electrons. Let us designate by A the dimensionless ratio of the 
energy transfered by the fast electron to the secondary electron, i.e., Ei — E[{= E 2 — m), 
which was originally at rest, to the kinetic energy of the fast electron Ei — mc^ = Ei — m. 
With this notation we can write the following formula for the differential cross-section of 


the electron-electron scattering (see, e.g., B) 


dcTp 


27ra(^ao 


dA nf(/5-l) A2(l-A)2 


{1 - [3 - (^) j A(1 - A) + - A)'^} (24) 


P - 1\2 

J 


where /3 = -A is the beta factor of the fast electron. It is interesting to note that the 
differential cross-section of the electron-electron scattering is a relatively simple function 
of the parameters j? and A (1 — A). For small A one finds • In dense 

media the energy loss due to electron-electron scattering is in direct competition with the 
bremsstrahlung and non-elastic collisions of fast electron(s) with atoms/ions. 

The last process which we want to discuss here is the Compton scattering of photons by 
a free (fast) electron. Actual high-temperature plasmas contains a large number of photons 
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which can interact with free electrons. Snch a process leads to the emission of secondary 
photons. The energy of the secondary photon can be larger, or smaller than the energy of 
the original photon. In the first case we deal with the Compton scattering of photon by a 
fast electron. The differential cross-section of the Compton scattering is written in the form 
(see, e.g., Q) 

dac ^ 


•‘0^2 





, 1 X 

/Ki 

^2\1 

4 ( 

H-) 

4 ( 

H-j 

- — + 

— 




K2^ 

^K2 

kiE 


(25) 


du2 m?Ki 

where U 2 = hu 2 is the energy of the secondary photon, while the parameters of the incident 
and final electrons ni and K 2 eqnal the prodncts of the electron and photon 4-vectors, i.e., 
cind K 2 = — respectively. As follows from the energy conservation law(s) 
for the Compton scattering, the same parameters Ki and K 2 can be written in the forms: 

^cid K 2 = — The freqnency of the secondary photon can de evalnated 

from the following relation for the fonr-vectors: pi + ki = P 2 + k 2 (energy conservation 
law). Indeed, calcnlating the both sides of the eqnality (pi -|- ki)"^ = {p 2 + ^ 2 )^ with the 
nse of the conditions Pi=P 2 = and kf = k 2 = 0 one finds the following identity 

Piki = Pik 2 + kik 2 , or in other words: 


U}lU}2 


a;i(l — vi cos 6^1) = a;2(l — vi cos6*2) H — ^ —(1 — vi cos6*) 

El 


(26) 


where vi and Ei are the velocity and energy of the incident electron, ui and 002 are the 
freqnencies of the incident and secondary photons, while 9i and 6*2 designate the angles 
between vectors pi and ki and vectors pi and k 2 , respectively. Analogonsly, the angle 9 
is the angle between ki and k 2 . More details abont snch calcnlations can be fonnd, e.g., 
in Section 3.7 of [^. Note that in high-temperatnre plasmas the Compton scattering of 
photons by fast electrons is one of the leading channels of energy loss. This process is 
always competing with the high-temperatnre bremsstrahlnng. 

A. Bremsstrahlnng from high-temperatnre plasma 


Bremsstrahlnng from hot plasma with temperatnres T ^ 7 - 15 keV is of great interest 
for actnal applications mainly related to nnclear fnsion. In general, snch a plasma can be 
considered as a system at local thermal eqnilibrinm. Therefore, we can investigate this by 
nsing a set of additional eqnations which follow from conditions of thermal eqnilibrinm and 


12 











allow us to obtain a number of relations between different properties of radiating plasma. In 
this case we can introduce separate temperatures for each plasma component, i.e. for ions T,, 
electrons T^. and radiation T^. To simplify the description of bremsstrahlung and/or Compton 
scattering below we consider a two component plasma which contains only electrons and 
radiation with temperatures Te and T^. It is clear that a hot plasma of light elements Q < 3 
with temperatures T 7 -15 keV contains only electrons and bare atomic nuclei. For such a 
plasma only bremsstrahlung and/or Compton scattering of electrons are important channels 
of energy loss. Without an extensive discussion of energy balance in fusion-related plasmas 
we just present the explicit formula for the coefficient A^r which determines the energy 
transfer rate (or temperature transfer rate) betweeen electrons and radiation in a high- 
temperature plasma which is assumed to be at local thermal equilibrium. If Tg and are 
the radiation and electron temperatures of the plasma, then we can write Aer = Cye{.^b + ^c)i 


where 0^^ = -- Tf 


15| . and the overall bremsstrahlung rate is 


= 8.510768649 • 


ke^N\ pZ 


he (rrigfcTg)2 




(27) 


where Na = 6.02214179-10^^ is the Avogadro number, a is the Stephan-Boltzmann constant, 
c is the speed of light in vacuum, rUg and e are the electron mass and electric charge of the 
electron, respectively, while k is the Boltzmann constant and G{x) is the following function 


G{x) = 


X 


_ 7+°° de/(0{l-expK(^-l)]} 

lio l-exp(-|) 


(28) 


where /(^) is the function 

/(O = exp(0 [ hi(^j/ + 1 + exp{-^y)dy 

^ 0 


(29) 


The formula, Eq. (|27)l . determines the ‘pure’ bremsstrahlung rate for Tg < Tg and inverse 


bremsstrahlung rate, if Tg < T^. The Compton scattering rate uq is 


m 


Uc = 134.0412866 • 


eV 


[rrieC 




kTt 


(30) 


Note that the product of the Boltzmann constant and Avogadro number equals the universal 
gas constant i?, i.e. kNA = R- 

These formulas can be used to determine a heat propagation in the light-element plasma 
from an initial (hot) area into new (cold) areas of the plasma. In general, if energy losses 


13 










during such a heat propagation into cold areas are over-compensated by the energy release 
from thermonuclear reactions, then this plasma ignites, or burns-up. In actual application 
the high-temperature burn-up of deuterium plasmas which contain signihcant amounts of 
deuterium and tritium nuclei and which can also be mixed with some amounts of the ^He 
and ®Li nuclei, is of great interest. Moreover, it is clear that the burn-up in high-temperature 
plasma will be optimal if the nuclear fusion reactions start in one very small spatial area (or 
‘point’) and propagate from this hot point to other areas occupied by the cold thermonuclear 
fuel. 


B. Bremsstrahlung role in the high-temperature burn-up of the deuterium plasma 


Let us assume that inside of an inhnite thermonuclear fuel with density po, part of the 
fuel in volume V is instantaneously heated to a high temperature T. If the values of V and 
T exceed some critical values Vc and then thermonuclear burning begins in the volume Vc 
and such a burning wave can propagate to the rest of the fuel [l6| - [l8|. This corresponds 
to the burn-up of thermonuclear fuel from a central hot point. In spherical geometry this 


process is governed by the following equation 

3 dvf 


16| 


dt 


-C 




Q{x,p,T) 

cV 

^ y m.n.T 


(31) 


where C is the specihc heat per unit mass of the fuel, r/(f) is the radius of the hot zone (also 
called the combustion zone), x = PqTj is the burn-up parameter and dx = podvf = poVmaxdt, 
where Vmax is the speed of burn wave propagation (at givien po and T), or velocity of the hot 
zone expansion. The notation Q{x,p,T) in Eq. fldip stands for the energy release function 
which depends upon chemical composition of the fuel, burn-up parameter x and density p. 
In actual high-temperature plasmas the burning wave can propagate in a few different ways, 
e.g., by compressing initially cold fuel by a very strong shock wave (or by a consequence of 
such waves), or by high-temperature heat conduction. The hrst way corresponds to the high- 
temperature detonation 17|], while the second way represents a high-temperature electron 
heat wave. In some plasmas at certain conditions a few other ways of heat propagation are 
also possible, e.g., by fast a—particles (or other electrically charged fast particles) formed 
in nuclear fusion reactions and/or fast neutrons which are often formed in such reactions. 
Below, we consider only actual plasmas in which propagation of the hot zone proceeeds by 
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high-temperature detonation and/or high-temperature thermal conductivity. 

The crucial question about any possible thermonuclear fuel is the explicit form of the 
burn-up function Q{x,p,T) in E g .03II) . For instance, for equimolar deuterium-tritium mix¬ 
ture this function takes the form 1 ^ 

Q{x,p,T) = (l + 4fc„)-2.17-10" 

31VT 


1 + 0.232r°-^5 , 20 , 

— 2 — , • exp(- 2 - 

TsVi + g.d- io-5r3-25 Ts 


(32) 


1 + i.i^r-i-75 

0.048X jg Qf energy which neutrons are leaving in a sphecial 


where K = i+o.o 48 x 

ball of radius r located inside of the equimolar DT-plasma with density p. The hrst term 
in Eq. (l3^ is the product of the fusion reaction cross-section and velocity of the colliding 
nuclei {av) averaged over Boltzmann energy distribution. For the (d, f)—nuclear reaction 
this (hrst) term is very large, since the energy realease from this reaction is huge and it has 
a very large resonance a.t T ^ 107 keV. For a pure deuterium plasma, the analogous term is 
substantially smaller (59 - 78 times smaller depending on temperature). The explicit formula 
of the hrst term in Eq. fl52D for the (d, d)—reaction can be found, e.g., in 19|. Therefore, for 
pure deuterium plasmas the second term in Eq. fl3^ becomes important. Breihy, the second 
term represents energy loss from the reaction (hot) zone due to hux of high-temperature 
radiation. This term is often called the bremsstrahlung (negative) contribution. In earlier 
astrophysical studies performed 70 - 75 years ago for hydrogen plasmas this term was written 
in the form —31\/T. However, such a form for the bremsstrahlung term in Eq. fl32|) leads to 
a uniform conclusion that ignition of the deuterium plasmas in hnite volumes with spatial 
radius less than 1 meter is not possible. Formally, thermonuclear ignition is possible when 
spatial radius of the ball of deuterium plasma exceeds 800 meters, but it is absolutely non- 
realistic in applications. After extensive reasearch of deuterium plasmas it became clear 
that all radiation quanta emitted from the hot zone due to high-temperature bremsstrahlung 
cannot reach the bondaries of the hot zone. This means that some part of these quanta is 
absorbed in this hot zone. In general, the number of absorbed quanta is directly proportional 
to the density p of the plasma and burn-up parameter x. The factor -|-1.1 ■ yTp- 
in the second term of Eq.([32]) represents the bremsstrahlung screening. Note that this factor 
contains two control parameters (density and linear size). Variation of these parameters and 
hrst of all density p allows one to produce ignition of pure deuterium and other thermonuclear 
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plasmas. This explains a crucial role of high-temperature bremsstrahlung for ignition of 
various deuterium-containing plasmas. 


IV. CONCLUSION 

We developed an approach which can be used to evaluate contributions of the electron- 
electron correlations into bremsstrahlung cross-section and corresponding energy loss for 
few-electron ions and atoms. The crucial part of accurate atomic calculatons is to determine 
the form-factors F{q) of the incident few-electron ions/atoms, i.e. the Fourier-transform of 
the one-electron density distribution. In the central held approximation (which is very 
good for these problems) actual computations are reduced to the accurate computation of a 
few one-dimensional integrals. Such computations have been performed for all two-electron 
ions/atoms (He - Ni^®’*') yielding the closed analytical expression for all matrix elements 
which are needed to determine atomic form-factors of all two-electron atoms and ions to 
very high accuracy within a central held approximation in which the electron distribution 
density is considered to be a spherically symmetric. There is a possibility to generalize our 
approach to three- and four-electron atoms/ions and even to the case of atoms/ions with 
arbitrary number of bound electrons. We also discuss the energy loss due to bremsstrahlung 
from high-temperature plasmas which contain multi-charged ions and free electrons. It is 
shown that screening of bremsstrahlung in high-temperature plasma plays a central role in 
the burn-up of high-temperature deuterium plasmas. 
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TABLE I: Convergence of the atomic form-factor in a.u. of the He-atom (the ground 1^S'—states) 
upon the total number of basis functions used. Atomic nucleus is assumed to be infinitely heavy. 


q 

F{q) {N = 2000) 

F{q) {N = 2500) 

F{q) {N = 3000) 

F{q) {N = 3500) 

F{q) (N = 4000) 

0.1 

1.99602833554455 

1.99602833554456 

1.99602833554457 

1.99602833554456 

1.99602833554456 

0.2 

1.96472444168180 

1.96472444168181 

1.96472444168181 

1.96472444168180 

1.96472444168180 

0.3 

1.93800126413211 

1.93800126413213 

1.93800126413214 

1.93800126413212 

1.93800126413213 

0.4 

1.90452526569111 

1.90452526569110 

1.90452526569111 

1.90452526569110 

1.90452526569110 

2.5 

0.80608623318909 

0.80608623318908 

0.80608623318909 

0.80608623318908 

0.80608623318908 

2.7 

0.72156548850375 

0.72156548850375 

0.72156548850375 

0.72156548850376 

0.72156548850375 

2.8 

0.68250273317700 

0.68250273317701 

0.68250273317702 

0.68250273317701 

0.68250273317701 

5.7 

0.01486064557919 

0.01486064557918 

0.01486064557919 

0.01486064557918 

0.01486064557918 

5.8 

0.01417307437410 

0.01417307437411 

0.01417307437411 

0.01417307437412 

0.01417307437411 

6.0 

0.01290641070551 

0.01290641070550 

0.01290641070551 

0.01290641070550 

0.01290641070550 
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TABLE II: Atomic form-factors in a.u. for some two-electron ions (ground 1^5—states). Atomic 
nuclei are assumed to be infinitely heavy. 


q 

F ( q ) ( H -) 

F { q ) ( He ) 

F { q ) ( Ne »+) 

F { q ) ( Cal »+) 

F { q ) ( Ni 26+) 

0.0 

2.0 

2.0 

2.0 

2.0 

2.0 

0.1 

1.96132945057312 

1.99602833554456 

1.999891483947858 

1.999973970809488 

1.999986873353467 

0.2 

1.85640934479775 

1.96472444168180 

1.999565989597304 

1.999895886309433 

1.999947494193332 

0.3 

1.71043482420304 

1.93800126413213 

1.999023678306149 

1.999765755713475 

1.999881864857887 

0.4 

1.54833795764938 

1.90452526569110 

1.998264818804848 

1.999583594374983 

1.999789989243940 

0.5 

1.38717714574940 

1.86490379021588 

1.997289786898069 

1.999349423783013 

1.999671872806301 

0.6 

1.23602867177643 

1.81982405765591 

1.996099065047703 

1.999349423783013 

1.999527522557058 

0.7 

1.09860493965672 

1.77002688319127 

1.994693241838125 

1.999063271556655 

1.999356947064652 

0.8 

0.97567464681582 

1.71628065247207 

1.993073011324763 

1.998335163282062 

1.999160156452743 

1.0 

0.77020482817221 

1.65935700342965 

1.989192627248489 

1.997399612788159 

1.998687978132901 

2.0 

0.25140513353668 

1.05558832506606 

1.957298839612961 

1.989629032969506 

1.994759690186353 

2.1 

0.22631640202285 

1.00140920768267 

1.953000585093706 

1.988570603016670 

1.994223728568467 

2.2 

0.20400578198386 

0.94931116662983 

1.948508086746243 

1.987461456215148 

1.993661855738936 

2.3 

0.18414561608885 

0.89936918205287 

1.943823468832696 

1.986301721988844 

1.993074104874982 

2.4 

0.16644818561710 

0.85162372319241 

1.938948936453435 

1.985091535520344 

1.992460510663290 

2.5 

0.15066077068448 

0.80608623318908 

1.933886773267480 

1.983831037715909 

1.991821109295450 

2.6 

0.13656133733834 

0.76274416108052 

1.928639339147500 

1.982520375169048 

1.991155938463219 

2.7 

0.12395476256023 

0.72156548850375 

1.923209067773552 

1.981159700122672 

1.990465037353574 

2.8 

0.11266952478171 

0.68250273317701 

1.917598464169758 

1.979749170429865 

1.989748446643585 

2.9 

0.10255480002725 

0.64549643651892 

1.911810102188170 

1.978288949513273 

1.989006208495098 

3.0 

0.093477913249790 

0.61047816022939 

1.905846621944145 

1.976779206323141 

1.988238366549214 

5.0 

0.0019119691603341 

0.02092499944253 

1.754339984527224 

1.936491784524877 

1.967584769495008 

5.1 

0.0017868988736314 

0.01990385784924 

1.745397409937382 

1.933990377128254 

1.966291880377667 

5.2 

0.0016716011208328 

0.01893974669421 

1.736349044142877 

1.931444462097725 

1.964974708365103 

5.3 

0.0015652000117109 

0.01802919975723 

1.727198296048866 

1.928854321129029 

1.963633329186683 

5.4 

0.0014669071827418 

0.01716896720854 

1.717948575709361 

1.926220240087431 

1.962267819865867 

5.5 

0.0013760122531516 

0.01635600319142 

1.708603291878175 

1.923542508939155 

1.960878258710472 

5.6 

0.0012918743968055 

0.01558745381000 

1.699165849608033 

1.920821421682095 

1.959464725302800 

5.7 

0.0012139148921932 

0.01486064557918 

1.689639647900089 

1.918057276275817 

1.958027300489612 

5.8 

0.0011416105304558 

0.01417307437411 

1.680028077405945 

1.915250374570904 

1.956566066371970 

5.9 

0.0010744877767607 

0.01352239490186 

1.670334518184135 

1.912401022237641 

1.955081106294943 

6.0 

0.0010121175936764 

0.01290641070550 

1.660562337512848 

1.909509528694104 

1.953572504837169 

10.0 

0.0001516730439791 

0.002624513248630 

1.24042038231789 

1.763227052540513 

1.874953401633500 

10.1 

0.0001459992755751 

0.002536906728319 

1.22994632121624 

1.758910183840519 

1.872561989249619 

10.2 

0.0001405861697800 

0.002452820055244 

1.21950396802771 

1.754566548761648 

1.870151446564242 

10.3 

0.0001354196056549 

0.002372089736317 

1.20909496362771 

1.750196545744651 

1.867721899795742 

10.4 

0.0001304863393311 

0.002294560815708 

1.19872090064025 

1.745800573630203 

1.865273475839531 

10.5 

0.0001257739428235 

0.002220086384021 

1.1^1^8332376428 

1.741379031587553 

1.862806302254074 

10.6 

0.0001212707475630 

0.002148527118019 

1.17808373013037 

1.736932319043730 

1.860320507246883 

10 7 

0 nom 1 f^QR.^7Q99FinFi 

n nn9n7Q7Kn84SS7s 

1 

1 

1 8Fi7S1fi91Qf^Rn4Ql 
















TABLE III: The electron-nuclear {r"^) expectation values (in a.u.) for a number of two-electron 
ions with infinitely heavy atomic nuclei. 


ion/atom 

H- 

He 

Ne®+ 

Ca33+ 

Ni2®+ 

i^eN) 

(rlN) 

{rl%) 

11.913699678051262 

645.1445424122194 

87266.1424069593 

22035718.695491 

8937970758.421 

1.19348299501893527 

3.9735649316629101 

26.28244697552571 

289.827674593716 

4797.7265440809 

3.2556160988739701-10-2 

2.690792104845699-10-3 

4.188297920816449-10-3 

1.05338028889520-10-3 

3.8992426943743-10-3 

7.8088339435014055-10-3 

1.5358087156342157-10-3 

5.6625739359113182-10-® 

3.364163448868867-10-’’ 

2.93640072611717-10-3 

3.93801344691995127-10-3 

3.8974100557929396-10-3 

7.222000786090979-10-’’ 

2.15474156619540-10-3 

9.5023197194145-10-’“ 
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